Operational Classification and Quantification of Multipartite Entangled States 
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We formalize and extend an operational multipartite entanglement measure introduced in T. R. 
Oliveira, G. Rigolin, and M. C de Oliveira, Phys. Rev. A 73, 010305(R) (2006) through the 
generalization of global entanglement (GE) [ D. A. Meyer and N. R. Wallach, J. Math. Phys. 
43, 4273 (2002)]. Contrarily to GE the main feature of this new measure lies in the fact that 
we study the mean linear entropy of all possible partitions of a multipartite system. This allows 
the construction of an operational multipartite entanglement measure which is able to distinguish 
among different multipartite entangled states that GE failed to discriminate. Furthermore, it is also 
maximum at the critical point of the Ising chain in a transverse magnetic field being thus able to 
detect a quantum phase transition. 
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I. INTRODUCTION 

Since Schrodinger's seminal paper 1} entanglement is 
recognized to be at the heart of Quantum Mechanics 
(QM). For a long time the study of entangled states was 
restricted to the conceptual foundations of QM 0, Q 
Since the last two decades, however, entanglement was 
also recognized as a physical resource which can be used 
to efficiently implement informational and computational 
tasks . The understanding of the qualitative and quan- 
titative aspects of entanglement, therefore, naturally be- 
came a fertile field of research. Nowadays, entanglement 
of bipartite states (a joint state of a quantum system 
partitioned in two subsystems A and B) is quite well 
understood. Good measures of entanglement for these 
systems are available, specially for qubits [f|. On the 
other hand entanglement of multipartite states (a joint 
state of a quantum system partitioned in more than 
two subsystems) cannot be understood through simple 
extensions of the tools and measures employed for bi- 
partite entangled states. Most of the tools available 
to study bipartite states (e.g. the Schmidt decomposi- 
tion 0) are in general not useful for multipartite states. 
Even a qualitative characterization of the many possi- 
ble multipartite entangled states (MES) is very com- 
plex since for a given A^-partitioned system there are 
many "kinds" of entanglement For example, let 

\^)n = |0i) • • • <8 \4> P } <S> \iP)n- p be a A^-partite state 
in which \4>i), 1 < i < p, is the ith subsystem state and 
\iP)m~ p is the state describing the other N — p subsys- 
tems. If \i/j)n~p is an entangled state then \*f?) is called 
a p-separable state 0. After discovering the value of p 
for a given multipartite state another complication shows 
up when we focus on \ip)N- p since its subsystems can be 
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entangled in several inequivalent ways. For example, in 
the case of three qubits there are two paradigmatic MES 
which cannot be converted to each other via local oper- 
ations and classical communication (LOCC)|3- For four 
qubits, nine different kinds of entanglement are possible, 
which cannot be converted to each other via LOCC 
Thus after considerable work we still lack a deep under- 
standing of MES and new tools must be developed in 
order to capture the essential features of genuine multi- 
partite entanglement (ME). 

Our aim in this paper is to shed new light on the way 
ME is characterized and quantified. We intend to do this 
by formalizing and extending an operational ME mea- 
sure introduced in Ref. [T^j. We emphasize that it is an 
operational measure in the sense that it is easily com- 
putable, even for a multipartite state composed of many 
subsystems. This new measure can be seen as an ex- 
tension of the global entanglement and we call it, from 
now on, the generalized global entanglement: E G . The 
generalized global entanglement has several interes ting 
features, two of which were already explored in Ref. |lOl: 
(i) in contrast to the global entanglement measure [ll| it 
can identify genuine MES and (ii) it is maximum at the 
critical point for the Ising chain in a transverse magnetic 
field. Another important aspect of Eq is the fact that 
it has an intuitive physical interpretation. We can relate 
it to the linear entropy of the pure state being studied as 
well as with the purities of the reduced n-party states ob- 
tained by tracing out the other N—n subsystems [T^ll3| . 

This paper is organized as follows. In Sec. [H] we for- 
mally define E^ and we extensively discuss a few im- 
portant properties satisfied by the generalized global en- 
tanglement. In Sec. IIIII we calculate Eq for the most 
representatives MES. This gives us a good intuition of 
the meaning of Eq and illustrates its usefulness. We 

(n) 

also compare E G with other measures available, high- 
lighting the main differences and the advantages and dis- 
advantages of each one. In the same section we use E G ' 
to quantify the ground state multipartite entanglement 
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of the one dimension (ID) Ising model in a transverse 
magnetic field. Finally, in Sec. I1VI we present our final 
remarks. 

II. GENERALIZED GLOBAL ENTANGLEMENT 

Global entanglement (GE) was firstly introduced in 
Ref. [llj to quantify the ME contained in a chain of N 
qubits. Latter it was demonstrated |14j to be equivalent 
to the mean linear entropy (LE) of all single qubits in the 
chain. This connection between GE and LE considerably 
simplified the calculation of GE and also extended it to 
systems of higher dimensions. An intuitive, though not 
so rigorous, way of understanding GE is to consider it 
as quantifying the mean entanglement between one sub- 
system with the rest of the subsystems. In this process 
we are dividing a system of N components into a sin- 
gle subsystem and the remaining N — 1 subsystems. We 
could, nevertheless, separate the system into two parti- 
tion blocks, one containing L subsystems and the other 
one N — L [Tfil Il7j . There are many different ways to 
construct a given "block". In Refs. [lijIlB] a block of L 
subsystems consisted of the first L successive subsystems: 
L = {Si, S2, S3, . . . , Sl}- But any other possible combi- 
nation of L subsystems could be employed to construct 
a block. We may have, for instance, a block formed by 
the first L odd subsystems: L = {Si, S3, S5, . . . , Szl-i}- 
It is legitimate to compute the LE of each one of these 
possible partitions. Roughly speaking this allows us to 
detect and quantify all possible 'types' of entanglement 
in a multipartite pure state. The generalized global en- 
tanglement (Eq^) is defined to take into account all of 
those possible partitions of a system composed of N sub- 
systems. Before we define E^ we highlight two of its 
main important qualities: (a) It is a relatively simple and 
operational measure. Since it is based on LE it can be 
easily evaluated and it is valid for any type of multipar- 
tite pure state (states belonging either to finite or infinite 

dimension Hilbert spaces); (b) Each class of Eq is re- 
lated to the mixedness/purity of all possible n-partite 



reduced density matrices out of a system composed of N 
subsystems, and thus it is not restricted to reduced den- 
sity matrices of only one subsystem as the original GE 
[llL Ha ] . This fact is helpful for the physical under- 
standing of E^ . 

(n) 

Following the definition of E G we move to the study 
of the general properties of this new measure relating it to 
the mixedness/purity of the various reduced density ma- 
trices of the system. After that we particularize to qubits 
focusing on the ability of the generalized global entan- 
glement to classify and quantify MES. We conclude this 
section by presenting a variety of examples, which clarify 

(n) 

the necessity of all the classes of E G , i. e. 71 = 1,2,3,..., 
to properly understand the many facets of MES. 

A. Formal Definition of the Measure 

Consider a system S which is partitioned into N sub- 
systems Si, 1 < i < N. Let \^f) 6 H be a quantum state 
describing S and H the Hilbert space of the whole sys- 
tem. Since we have N subsystems, TL = Hi ® ■ • -®TLn = 
m which Hi is the Hilbert space associated with 
Si . The density matrix of S is p = | ^)(^! | and we define 
the generalized global entanglement [TJj as, 

JV-l N-l N-l 

4 k) (p) = E E ••• 

n— 1 j 1= l i 2 =ii+l 13=22+1 
N-l 

E G(n,h,i 2 ,...,in-i), (1) 

i n _i=i n _2+l 

where all the parameters are natural numbers, n < N, 
and 

N _! = (N-iy. 

n - x (N-n)\(n-l)\ 

is the definition of the binomial coefficient. Note that 
the summation is over all ikS, with the restriction that 
1 < ii < i 2 < ■ ■ ■ < 2,1-1 < N — 1. We also assume 
io = 0. The function G is given as, 



N- 



G(n,ii,i 2 , 



1 - 



TV- 



E 

3=1 



Tr 



(2) 



where Pj,j+i 1 ,j+i 2 ,...,j+i n - 1 is obtained by tracing out all 
the subsystems but Sa = {Sj, Sj^ , Sj+i 3 , . . ■ , Sj+i n _^} 
and d = min{dimS' J 4,dimS' J 4}. Here dimSU and dim 5^ 
are, respectively, the Hilbert space dimension of the sub- 
system Sa and of its complement Sa- I n resume the 
index n is for the number of subsystems in the A parti- 



tion and the indexes ii,i 2 , ■■■,i n ~i are the neighborhood 
addressing for each of the involved subsystems. 
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B. General Properties 

Eqs. Q and 10) are valid for any multipartite pure 
system, even systems described by continuous variables 
(Gaussian states for example). The key concept behind 
generalized global entanglement is the fact that it is 
based on the linear entropy, which is an entanglement 
monotone easily calculated for the vast majority of pure 
states. Thus, by its very definition, Eq and G inherit 
all the properties satisfied by LE, including the crux of 
all entanglement monotones: non-increase under LOCC. 

Another important concept of Eq^ and G is the in- 
troduction of classes of multipartite entanglement (ME) 
labeled by the index n. As we will see, they are all re- 
lated with the many ways a multipartite state can be en- 
tangled. Moreover, a genuine n-partite entangled state 
must have non-zero Eq and G's for all classes n. Here 
a genuine MES means a multipartite pure entangled sys- 
tem in which no pure state can be defined to anyone of 
its subsystems. There is only one pure state describ- 
ing the whole joint system. For three qubits, for in- 
stance, the states \GHZ) = (l/\/2)(|000) + | 111)) and 
\W) = (1/V3)(|001) + |010) + |100>) are genuine MES 
but |f) = (l/\/2)(|00) + |11))|0) is not. 

(n) 

Let us now explicitly show how the first classes of E G 
look like. This will clarify the physical meaning of the 
measure as well as the intuitive aspects which led us to ar- 
rive at the general and formal definitions given in Eqs. ifTJl 
and JSJ. 



1. First Class: n = 1 
When n = 1 Eqs. and J2J) are the same, 



4 1) ( P ) = G(1) = 



1 



3=1 



(3) 



and if we remember the definition of the linear entropy 
for the subsystem j [lj, [lj| , 

E L (p j ) = 1 ? ri [l-Tr(fi})], (4) 
then Eq. © can be written as 01 
1 N 

E G=„ 1 EE L (p j ) = (E L (p j )). (5) 

3=1 

In other words, Eq^ is simply the mean linear entropy of 
all the subsystems Sj . We should mention that for qubits 

(d = 2), Eq^ was shown to be exactly the Meyer and 
Wallach global entanglement [Tl|. 

The physical intuition behind the study of the mean 
linear entropies lies in the fact that the more a state is a 



genuine MES the more mixed their reduced density ma- 
trices should be. However, we should not limit ourselves 
to evaluating the reduced density matrices of single sub- 
systems Sj. We can take either two, or three, or n 
subsystems and calculate their reduced density matrices 
and also calculate their mean linear entropies. This is 
the reason of why we introduced the other classes of gen- 
eralized global entanglement. 



2. Second Class: n = 2 

For n = 2 Eqs. and © are not identical anymore, 
being, nevertheless, entanglement monotones: 



4 2) (p) 



N-l 



G(2,h) = 



N - 



h E 0(2, 



»i=l 



1 



N-n 



3=1 



(6) 



•(7) 



Now we deal with the reduced joint density matrix for 
subsystems Sj and Sj^ . The extra parameter i\ is 
introduced to take account of the many possible 'dis- 
tances' between the two subsystems. For nearest neigh- 
bors i\ = 1, next-nearest neighbors «2 = 2, and so forth. 

Noticing that the linear entropy of the subsystems Sj 
and Sj+ij by tracing out the rest of the other subsystems 
is given by 



E L (p j>j+il ) = — [l-Tr 



(8) 



then Eq. Q can be written as, 
This implies that Eq. © is simply given by 



(9) 



4 2) (p) 



N-l 



N 



L-J2( E L(Ps,j+h)) = ((EL(Pj, j+il ))), 



(10) 

where the double brackets represent the averaging over 
all possible G(2, 1 < i x < N - 1. 

Looking at Eqs. 10 and 11U|) we can easily interpret 
E ( e ] and G(2,ii). First, let us deal with G(2,h). We 
assume that all the subsystems are organized in a linear 
chain. (This assumption simplifies the discussion in what 
follows.) If we remember that l<«i<JV— 1, where N is 
the number of subsystems, Eq. JJjJ tells us that G(2,ii) 
is nothing but the mean linear entropy of two subsystems 
with the rest of the other subsystems conditioned on that 
these two subsystems are i\ lattice sites apart. 

For concreteness, let us explicitly write all the possible 
G(2, i\) for a linear chain of five subsystems. Since N = 5 
we have 1 < i± < 4, which gives four G's pictorially 
represented in Fig. ^ 
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(1) G(2,l), which is the mean linear entropy (LE) of 
the following pairs of subsystems with the rest of 
the chain: {(Si, S 2 ), (S 2 ,S 3 ), (S3, Si), (£4, £5)}; 

(2) G(2, 2), which is the mean LE of the following pairs 
of subsystems: {(Si, S3), (£2, £4), (S3, £5)}; 

(3) G(2, 3), which is the mean LE of the following pairs 
of subsystems: {(Si, S4), (S 2 , £5)}; 

(4) G(2, 4), which is the mean LE of the following pairs 
of subsystems: {(£1, S5)}. 



(!)••••• 

Sj + S 2 S 3 + S 4 



(2) O O O O O 



Sj + s 3 s 3 + s 5 



(3) O O O 



Si + s 4 



(4) 



Sj + S, 

Figure 1: All combinations of two elements out of five. 

Finally Eq. shows that E£> is the mean linear 
entropy of two subsystems with the rest of the chain ir- 
respective of the distance between the two subsystems, 
i.e., it is the averaged summation of all the (l)-(4) kinds 
of G(2,h), 1 < h < 4. 

3. Third Class: n = 3 
By setting n = 3 Eqs. and J2J become 

N-l N—l 



4 3) w = 



and 



(N — l)(N -2) ^ ^ 

v n ' ii=l i 2 =ii+l 



£ E G(n,h,i 2 ), (11) 



G(3,ii,i 2 ) = 



d- 1 



N — Z2 



3+12 1 



(12) 

Eq. I|12|l deals with reduced density matrices of three sub- 
systems: Sj, £j+i 15 and £/+i 2 . Therefore, G(3, 11,^2) is 



the mean linear entropy of all three subsystems with the 
rest of the chain conditioned to that Sj+i 1 and £j+j 2 are, 
respectively, i\ and i 2 lattice sites apart from Sj. Taking 
the mean of all possible G(3,ii,i 2 ) we obtain Eq. lllfl . 
This is equivalent to averaging over all linear entropies 
of three subsystems irrespective of their distances. Al- 
though we do not explicitly write them here, similar ex- 
pressions as those given by Eqs. 10 and I|1U|I can be ob- 
tained for this class. 

Again, as we did for the second class, it is explanatory 
to analyze in details the N — 5 case. Now 1 < ii < i 2 < 
4. This time we have six G's (See Fig.|2J: 

(1) G(3, 1,2), which is the mean linear entropy (LE) of 
the following triples of subsystems with the rest of 
the chain: {(Si,S 2 ,S 3 ), {S 2 ,S 3 , £4), (£3, £4, £5)}; 

(2) G(3,l,3), which is the mean LE of the following 
triples of subsystems: {(£1, £2, £4), (£2, £3, £5)}; 

(3) G(3, 1,4), which is the mean LE of the following 
triples of subsystems: {(£1, £2, £5)}; 

(4) G(3,2,3), which is the mean LE of the following 
triples of subsystems: {(£1, £3, £4), (£2, £4, £5)}; 

(5) G(3,2,4), which is the mean LE of the following 
triples of subsystems: {(£1, £3, £5)}; 

(6) G(3, 3,4), which is the mean LE of the following 
triples of subsystems: {(£1, £4, £5)}. 



s 2 + s 3 + s 4 



s 2 + s 4 + s 5 



r 



I J J 

(i) • • • • • <«> • • • 1 

t t ft t t t t 1 

s, + s 2 + s 3 s 3 + s 4 + s 5 s, + s 3 + s 4 

s 2 + s 3 + s s 



1 ] 



< r 



1 



(5) 



o o 



t L 



Sj + s 2 + s 4 



Sj + S 3 + S s 



(3) o O O O O («) o o o o o 

t t t t 1 f 

s, + s, + s, s, + s 4 + s, 

Figure 2: All combinations of three elements out of five. 



4- Higher Classes: n > 4 

Remembering that n < N, higher classes n of E < q\p) 
only make sense for systems such that N > n + 1 sub- 
systems. The higher a class n the greater the number 
of G's necessary for the computation of E@\p). This is 
a satisfactory property we should expect from a useful 
multipartite entanglement measure since as we increase 
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the number of partitions of a system we increase the way 
it may be entangled 0, IU • 

If we employ the definition of LE for n subsystems out 
of a total of N, 



El( P . 



j,...,j+t n -± l 



d-1 



1 - Tr 



we can write Eqs. and (JIJ respectively as 



G(n,h, . 



■ .,ln-l) 

E { a\ P ) 



( E L (pjj+ii,.. 

({E L (p jij+il! 



A.-J), 



(13) 

(14) 
(15) 



In Eq. 114|) the single pair of brackets ( ) represents the 
averaging over all possible configurations of n subsystems 
in which subsystem Sj+i k is ik lattice sites apart from Sj. 
Here 1 < k < n — 1. Finally, the double brackets (( )) 
is the average of the linear entropy of n subsystems over 
all possible combinations (distances) in which they can 
be arranged. 

We should mention at this point that Eq and G 

are more general than the block entanglement (E^ ) as 
presented in Refs. [l6L UtI Il9j. By block entanglement 
it is understood that we divide a set of N subsystems 
{Si, S 2 , ■ ■ ■ , S N } in two blocks, A n = {Si, S 2 , . . . , S„} 
and Bpf- n — {S n +i, S n +2) . . . , Spf}, and calculate the 
linear or von Neumann entropy between blocks A n and 
-Bjv-n- In the language of generalized global entangle- 
ment, block entanglement for a translational symmetric 
state is simply 



E 



(n) 



G(n,ii = l,i 2 = 1,. . .,in-i = 1), 



which is only one of the many G"s we can define. The 
main difference between these two measures lies in the 
fact that we allow all possible combinations of n subsys- 
tems out of N to represent a possible 'block'. Contrarily 
to block entanglement, here there exists no restriction 
onto the subsystems belonging to a given 'block' to be 
nearest neighbors. They lie anywhere in the system's 
domain. 



C. Particular Properties for Qubits 

Although Eqs. |QJ and © are defined for Hilbert 
spaces of arbitrary dimensions we now focus on some 
properties of Eq and G for qubits. There are two main 
reasons for studying qubits in detail. Firstly, they are rec- 
ognized as a key concept for quantum information theory 
and secondly, the simplest multipartite states are con- 
structed employing qubits. 

Let p = \^)(^f\ be the density matrix of a N qubit 
system and pj = Try(p) the reduced density matrix of 
subsystem Sj, which is obtained by tracing out all sub- 
systems but Sj . A general one qubit density matrix can 
be written as 



Pj = Trj(p) 



(16) 



where the coefficients are given by 

tf = Tr(<T? Pj )=(*\o?\*). 



(17) 



Here u" is the Pauli matrix acting on the site j, a = 
0, x, y, z, where a is the identity matrix of dimension 
two, and p" is real. Since pj is normalized po = 1. Using 
Eqs. and dJ) we obtain 

Tr ( Pj 2 ) = i (1 + (<r*) 2 + (aj) 2 + ( CT |> 2 ) . (18) 

This last result implies that Eq. @ can be written as 

1 N 

4 1} = i-^E(K) 2 + K) 2 + K) 2 )- (w) 



One interesting situation occurs when we have transla- 
tional invariant states p. (The Ising model ground state 
for example.) In this scenario (erf) = (cr") for any i and 
j. Therefore, Eq. I|19(l becomes 



23g> = 1 - (aj) 2 - (a}) 2 - <af> 2 , 



(20) 



which is related to the total magnetization M of the sys- 



1- 



tern, |M| 2 = NiW+W+W), by E, 

By tracing out all subsystems but Si and Sj we obtain 
the two qubit reduced density matrix 



where 



af3 



Tr (of of pa) = <*kf°?l*> 



(21) 



(22) 



Eq. (|21|) is the most general way to represent a two-qubit 
state and together with Eq. (1221) imply that 



4 

a,/3 



(23) 



Remark that in Eq. <\2',i\) the trace of pfj is the sum of all 
one and two-point correlation functions. Moreover, since 
E G and G{2,h) depend on Eq. we find in these 

entanglement measures both diagonal and off-diagonal 
correlation functions. 

Again it is instructive to study translational symmetric 

a/3 _ /3a 

°ij — Pij 
assumption in Eq. (Q we get 

2 



states in which p^ = p 1 :- for any a and (3. Using this 



G(2,ii 



1-3 [«> 2 + K) 2 + K) 2 + (^ y 



3 u j+iii 



33 

JV -V \2 



+ (« + u) 2 /2 + K^ 1 > 2 /2] 



(24) 



Note that the previous formula is not valid for N < 3. 
For N = 2 only E^ is defined and for N = 3 we have 
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d = minjdim Sa, dim Sa} — 2 and not d = 4, the value 
of d for all N > 4. Now if we compare G(2, ii) with the 
concurrence (a bipartite entanglement monotone), as we 
do for the Ising model in Sec. IIIIBI we will note that 
while the concurrence does not depend on any one-point 
and on any off-diagonal two-point correlation function 
[HUH G(2,h) does. 



D. Why Do We Need Higher Classes? 

The simple fact that different types of entanglement 
appear as we increase the number of qubits (or equiv- 
alently the number of subsystems) @, @ indicates that 
the various classes here introduced may be useful to clas- 
sify and quantify the many facets of ME. For example, 
the first class does not suffice to unequivocally quan- 
tify MES. Although it is maximal for Greenberger-Horne- 
Zeilinger (GHZ) states [2^ it is also maximal for a state 
which is not a MES, as we now demonstrate. Let us 
compute E(V for three paradigmatic multipartite states. 
The first one is the GHZ state: 



\GHZ N ) = ±(\0)® N + \1)® N ) 
v2 



(25) 



where 10)®^ and \1)® N represent, respectively, N tensor 
products of the states |0) and |1). The GHZ state is a 
genuine MES since by measuring only one of the qubits 
in the standard basis we know exactly the results of the 
other TV — 1 qubits. Furthermore, tracing out any one of 
the qubits we obtain a separable state. A direct calcula- 
tion gives E { e ] (GHZ N ) = 1. 

The second state we shall analyze is given by a ten- 
sor product of N/2 Einstein-Podolsky- Rosen (EPR) Bell 
states 12] : 



\EPR 



N 



|4> H 



(26) 



where |<J>+) = (l/\/2)(|00) + |11)). For definiteness, we 
chose one specific Bell state. However, the results here 
derived are quite general and valid for any N/2 tensor 
products of Bell states. This state is obviously not a 
genuine MES. Only the pairs of qubits (2j — 1, 2j), where 
j = 1,2, ...,7V, are entangled. Nevertheless, we again 
obtain Eq\ePRn) = 1- This last result illustrates that 
Eq^ being maximal is not a sufficient condition to detect 
genuine MES. Note that E$ for both the GHZ N and 
EPRn states are independent of the number of qubits 
N in the chain. 

The last state we consider is the W state Q- It is 
defined as, 



\W N ) = 



1 N 

—= V|000---l 7 ---000). 

.7 = 1 



(27) 



As shown in Ref. 
that 
(N 



E { g\w n ) 



Note 



4(A- 1)/A 2 . 

j g depends on N and at the thermodynamic limit 

► oo) we have E^\Wn) = 0. For three qubits, 
the W state was shown 7) to be a genuine MES not 
convertible via LOCC to a GHZ state. 



(2) 

The computation of E G ' and G(2, 1) give different val- 
ues for each of those states. Remark that for N — 2 the 
previous functions are not defined and that for N = 3 
E ( G } = G(2,l) = 1. Table ID shows E ( G ] and (7(2, 1) for 
the states GHZn,EPRn, and Wn- We should mention 

(21 

that due to translational symmetry, G(2, 1) and E G are 
identical for the states GHZm and Wn- It is interesting 



Table I: The third and fourth columns give G(2, 1) and Eq ' 
for the three states listed in the first column when N > 3. 
The second column gives E(P for all N. Contrary to Eq\ 
we see that G(2, 1) and Eq distinguish the three states from 
each other. 



E. 



Ci) 



G(2,l) 



E, 



(2) 



GHZm 
EPR n 

W n 



4(jV-l) 
N 2 



N-2 
2(JV-1) 

16(iV-2) 
37V 2 



2 
3 

(2N-l)(N-2) 
2(iV-l) 2 



16(AT-2) 
37V 2 



to note that depending on the value of N, the states are 
differently classified through G(2, 1). Fig. illustrates 
the behavior of G(2, 1) for those three paradigmatic state 



as we vary N. A similar behavior is observed for E t 



(2) 



— GHZ N 

-A-W., 



_o — o — o — o — o- 



The state 1 000 • • • 1 3 ■ • ■ ■ 000) represents a N qubit state 

in which the j'-th qubit is |1) and all the others are |0). the definition of Eq' . For the EPRn state, G(2,l) 



4 6 8 10 12 14 16 18 20 

N 

Figure 3: (Color online) Here we show G(2, 1) as a function 
of the number of qubits N for the states GHZn, EPRn and 
Wn- Note that only when N = 4 we have two states with 
the same entanglement. Furthermore, for 4 < N < 8, Wn 
is more entangled than EPRn- This ordering is changed for 
TV > 9. 

(Fig. 0J . In this case, however, EPRn is the most en- 
tangled state for long chains. The reason for this lies in 

(2) 
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1 for any I > 2. Therefore, since Eq" 1 is obtained av- Consider now the state defined as [2c 
eraging over all G(2, 1), for long chains G(2, 1) does not 



contribute much and E, 



(2) 
G 



1. 




Figure 4: (Color online) Here we show E^ as a function 
of N. Again, only when N = 4 we have two states with the 
same entanglement. Moreover, for N > 4, EPRm is the most 
entangled state. 



We also calculated the values of ijff , E^ , and G(2, 1) 
at the thermodynamic limit. See Tab. [Hi Thus even at 



l<?l) = 2(l°)ll )2l >3|0)4+|0) 1 |l) 2 |0) 3 |l)4 
+|l)l|0) 2 |l) 3 |0)4+|l)l|l) 2 |l) 3 |l)4) 
= ^(|0)l|0) 3 |0) 2 |0)4+|0) 1 |0) 3 |l)2|l)4 

+ |l)i|l) 3 |0) 2 |0) 4 + |l) 1 |l) 3 |l) 2 |l) 4 ) 

= |* + )l 3 | < i ) + )24, (29) 

which is also a pair of EPR states. This time, however, 
subsystem Si is entangled with S3 and subsystem 5 2 
is entangled with S4 (See Fig. |SJ. Although different 



Alice 



Bob 



EPR, 



Alice Bob 



Figure 5: Pictorial representations of the states (a) EPR2 
and (b) g±. 



Table II: Eq\ G(2, 1), and E^ at the thermody namic limit. 



N -> 00 


^G 


G(2,l) 


g(2) 
G 


GH Zn 


1 


2/3 


2/3 


EPR n 


1 


1/2 


1 


Wn 












pairs of subsystems are entangled in these two differ- 
ent states, their amount of entanglement is the same: 
there are two EPR states in both cases. This fact is 
captured by the entanglement measures here introduced, 
i.e. E ( q\ePR 2 ) = E ( £\gi). The block entanglement, 
nevertheless, does not always give the same value for the 
two states above (see Tab. Illlll . This example illustrates 



the thermodynamic limit E G and G(2, 1) distinguish the 

three states. The ordering of the states, nevertheless, is 

(2) 

different. Again this is related to the definition of E G 
and is due to the contribution of G(2,l), I > 2, in the 
calculation of E { G ] {EPR N ). 

Besides a measure of multipartite entanglement be- 
ing able to distinguish different kinds of states it should 
not differentiate states that essentially contain the same 
amount of entanglement. For example, let us consider 
the following state, 

\EPR 2 ) = |$ + )i 2 |$+) 34 

= ^(|0)i|0) 2 |0) 3 |0) 4 + |0) 1 |0) 2 |l) 3 |l) 4 
+ |l)i|l) 2 |0) 3 |0) 4 + |l) 1 |l) 2 |l) 3 |l) 4 ).(28) 

This state describes a pair of EPR states where subsys- 
tem Si is entangled with 5 2 and S3 is entangled with S4. 



Table III: Comparison between Eq , G(2, 1), and Eg 





4 15 


p(2) 
G 


G(2,l) 




p (2) 
B 


EPR2 
91 


1 
1 


7/9 
7/9 


1/3 
1/3 


1 
1 




1 



that the block entanglement, as its name suggests, quan- 
tifies only the entanglement of partition A (sites 1 and 2) 
with partition B (sites 3 and 4). The generalized global 

entanglement E^, however, quantifies the amount of en- 
tanglement of a state independently on the way it is dis- 
tributed among the subsystems. We can go further and 
show the importance of using higher classes Eq to cor- 
rectly quantify the entanglement of a multipartite state 
no matter how the entanglement is distributed among 
the subsystems. For example, consider the state 

\GHZff) = \GHZ N )® M , (30) 

where the integer M > 1 represents how many tensor 
products of GHZn we have. Restricting ourselves to 



8 



N = 3 and M = 2 we get, 



1 



\GHZ() = -p(|000) + |lll))®-=(|000) + |lll)) 



V2 



V2 



-(1000000) + 1000111) + |111000) 

+1111111)). 



(31) 



Here, subsystems Si, S 2 , and S3 form a genuine MES and 
S4, S5, and Se another one. For this state E { g\GHZ$) = 
0. If we interchange the second qubit (S 2 ) with the fifth 
one (55) wc obtain the following state: 

\ZHG\) = i(|000000) + 1010101) + 1101010) 

+ 1111111)). (32) 

Now subsystems Si, S3, and S5 form a genuine MES 
and S 2 , 5*4, and Sq another one (See Fig. EJ). Those 



Alice 

1 S, 



(a) O S 2 



(bl 




GHZ, 



Figure 6: (Color online) Pictorial representations of the states 
(a) GHZl and (b) ZHG%. 

two states have the same amount of entanglement, i. e. 
two GHZ states. However, the computation of the block 
entanglement gives E { g\ZHG 2 3 ) = 6/7 ^ E ( g\GHZ 2 ). 
Had we employed the generalized global entanglement 
we would have obtained E { *\GHZl) = E%\ZHG%) in- 
stead. In general we have E^(GHZ^) ^ E { g\ZHG 2 n ) 

and E { q ] (GHZI) = E^ ] {ZHGl). Therefore, if we want 
to study the amount of entanglement of a multipartite 
state, independently on how it is distributed among the 

subsystems, we should employ Eq instead of E^\ since 
the later furnishes only the amount of entanglement be- 
tween a particular two block-partition in which the sys- 
tem can be divided. 



III. USEFULNESS OF THE GENERALIZED 
GLOBAL ENTANGLEMENT 

In this section we present two examples in which we 
explore the ability of E^' and the auxiliary measure 
G(n, ■ ■ • , in-i) to quantify multipartite entangle- 
ment. The first example deals with a finite chain of four 



qubits. We show that E G together with G(2, 1) allow 
us to correctly identify MES. Moreover, comparing the 
values of G(2,i\) for all the MES here presented we are 
led to a practical definition of what is a genuine MES. 
In the second example we investigate the entanglement 

properties of the Ising model ground state. We show 

(2) 

that Eq and G(2,ii) are maximal at the critical point 
and we analyze what correlation functions are responsible 
for this behavior of the generalized global entanglement. 
The results herein presented suggest that the long range 
correlations in the critical point for the Ising model are 
related to genuine MES. 



A. Finite Chains 

Let us now focus on the simplest non-trivial spin- 1/2 
chain, i. e. states with N = 4 qubits, by stud ying the 
entanglement pro perties of four genuine MES |24l l25| . 
The first one [24| is the famous four qubit GHZ state 

\GHZi) = |*i) = -L (10000) + |1111>) . (33) 
v2 

Qualitative and quantitative features of this state were 
already discussed in Sec. Ill Dl A direct calculation gives 



= 1;-®g ; (*i) = G(2,ii)($i) = 2/3, where 



r(2) 



'1 



1,2,3. The second state |24| is written as, 



|* 2 ) = ^=(V2|1111) 
+ 10001)). 



1 1000) + 10100) + 10010) 



(34) 



Calculating its first and second order generalized global 

entanglement we obtain Efr ($2) = l;i?cy ($2) = 
G(2,i 1 )($ 2 ) = 8/9. Note that as well as |$ x ) this 
state is a translational symmetric state. Moreover, 
G(2,ii)($ 2 ) > G(2,«i)($i). This last result will turn 
out to be very useful in constructing an operational def- 
inition of MES. The third state is given as, 



l*a) = -(11111) 



|1100) + |0010) + |0001)). (35) 



Since this state is not translational symmetric, G(2,ii) 
are not all equal. After a straightforward calculation we 

obtain Eg ($3) = 1;4 2) ( $ s) = 25/27; G(2,l)($ 3 ) = 
7/9;G(2,2) ($ 3 ) = G(2,3)($ 3 ) = 1. Again we should 
note that G(2,ii)($ 3 ) > G(2,ii)($i). 

These three states have in common a few remarkable 
properties j23| : (a) The local density operator describing 
each qubit is the maximally mixed state (l/2)i2, where J 2 
is the 2x2 identity matrix, thus explaining why E^ = 1 
for all of them, (b) The two- and three-qubits reduced 
operators do not have any fc-tangle j2||, k — 2,3. This 
emphasizes that they all are genuine MES, i. e. there is 
no pairwise or triplewise entanglement, (c) They cannot 
be transformed into one another by LOCC. 
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We shall consider a fourth state, 

\x) = ^(10000) -10011) -|0101) + |0110) 

+ |1001) + |1010) + |1100) + |1111)) , (36) 

recently introduced and extensively studied in Ref. |25j . 
The main feature of this state lies in its usefulness to 
teleport an arbitrary two-qubit state. Employing x this 
task can be accomplished either from subsystems S\ and 
52 to S3 and S4 or from S\ and S3 to S 2 and S4. The 
usual channel (two Bell states) used to teleport an ar- 
bitrary two-qubit state |23l |27| can teleport two qubits 
only from a specific location to another one: from Si 
and S2 to S3 and S4 for example. In addition state 
|x) has a hybrid behavior in the sense that it resem- 
bles both the GHZ and W states [25|. Tracing out 
any one of the qubits the remaining reduced density ma- 
trix cr has maximal entropy, a characteristic of the GHZ 
state. However, a has a non-zero negativity |28j | between 
one qubit and the other two [25[, a property of the W 
state. By calculating the generalized global entanglement 

we obtain ( X ) = lj^fr) = 23/27; G(2, 1) (x) = 
8/9; G(2, 2) ( X ) = 1; G(2, 3) (x) = 2/3. Again we see that 
for all h we have G(2+i)(x) > G(2+i)($i). 

We have grouped in Tab. II VI the entanglement calcu- 
lated for the previous four states. It is clear then that 



Table IV: Calculated values of Eq and G(2, i\) for the gen- 
uine MES shown in Sec. Ill I Al and for the EPR2 state. 





G 




G(2,l) 


G(2,2) 


G(2,3) 


EPR2 


1 


7/9 « 0.778 


1/3 


1 


1 


$1 


1 


2/3 « 0.667 


2/3 


2/3 


2/3 


$2 


1 


8/9 « 0.889 


8/9 


8/9 


8/9 




1 


25/27 « 0.926 


7/9 


1 


1 


X 


1 


23/27 « 0.852 


8/9 


1 


2/3 



Eq cannot be considered as the last word concerning 
the quantification and classification of MES. A glimpse 
of the first column in Tab. I1VI shows that all the five 
states listed have E G ^ = 1, even the EPR 2 state, an ob- 
vious non-genuine MES. Therefore, since E ( G } = 1 is not 
useful to classify different genuine MES or to correctly 
identify them we are compelled to go further and study 
the higher classes of the generalized global entanglement 

in order to achieve such a goal. Turning our attention to 

( 2) 

Eq ' we see that it is different for all the five states listed 

in Tab. IIVI implying that E G can distinguish among 

(2) 

the five states. According to E G the most entangled 
state is $3, which was shown to be a genuine MES [24(. 

Moreover, important clues for the understanding of 
what kind of entanglement is present in a given multi- 
partite state are also available in G(2, i\ = 1,2,3. 
Actually, these auxiliary entanglement measures give us 
a more detailed view of the types of entanglement a 

(2) 

state has than Eq since the latter is an average over 



all G(2+i). For example, if we relied only on E G to 
decide whether or not a state is a genuine MES we 

would arrive at a wrong answer. This point is clearly 

(2) 

demonstrated if we compare E£ J for the states EPR 2 
and $1 (GHZ4). Looking at Tab. El we see that 
E { G \EPR 2 ) > £^($1), where EPR 2 is not a ge nuine 

MES. The averaging process, as explained in Sec. Ill Dl 

(2) 

is responsible for this relatively high value of E G for 
the state EPR 2 . Remark that for translational symmet- 

ric states E G and G(2,ii) are equivalent to detect a 
genuine MES. However, if we analyze all the G(2,ii) 
terms we are able to detect a common characteristic 
shared only by the genuine MES: for x and all i\ we have 
G(2+i)($j,x) > G{2,i 1 ){GHZ i ) = 2/3. This suggests 
the following operational definition of a genuine MES: 

Definition 1 Let \*Sf) be a pure state describing four 
qubits. Z/G(l) = 1 and <?(2,ii)(tf) > G(2,i 1 )(GHZ i ) = 
2/3, i\ — 1,2,3, then |\&) is a genuine MES. 

Besides being practical, Definition^has a simple phys- 

(2) 

ical interpretation if we remember that E G and G(2, i\) 
are constructed in terms of the linear entropy of any two 
qubits with the rest of the chain. Noticing that the linear 
entropy is related to the purities of the two-qubit reduced 
density matrices, the definition above establishes an up- 
per bound for all the two-qubit purities of a MES. In 
other words, if all the two-qubit purities are below this 
upper bound the iV qubit state can be considered a gen- 
uine MES p9|. Furthermore, this upper bound was cho- 
sen to be that of the GHZ state, which is undoubtedly 
a genuine MES. 

Remark also that since G(2,ii) is a monotonically de- 
creasing function of the purities, an upper bound for the 
purities implies a lower bound for the value of G(2,ii) 
(cf. Definition 2]) . We can easily generalize this defini- 
tion to N qubits if we express it in terms of all n-qubit 
purities (n < N): 

Definition 2 A pure state of N qubits is a genuine 
MES if 

Tr(pl) < Tr(al) = 1/2, 
Tr(p* uj2 ) < Tr(al 2 )= 1/2, 

''■(I'l , > < M<2,..,„) = 1/2, 

where 

p,\j*,..jn = z*snra»(l*>(*l). 

= Tr V2 — ri {\GHZ N )(GHZ N \), 

and 

1 < h < n, 

1 < h < 32 < N, 
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l < h < h < • ■ • < jn < N. 



Note that as we increase the size of the chain we need 
to calculate more and more purities. Take for instance 
the state GHZ% given by Eq. A direct calculation 

gives Tt(pI) = 1/2 for I < ji < 6, Tr(p| 4 ) = 1/4, 
and Tr(p 3 2 J = 1/2 for all 1 < j 1 < j 2 < 6 but 
{jiiJi) — (3,4). Hence, if we restricted Definition |2 
just to the one- and two-qubits reduced density matrices 
we would erroneously conclude that GHZ^ is a genuine 
MES. Extending, however, the definition to all possible 
reduced density matrices we can detect that GHZ\ is not 
a genuine MES since Tr(p^ 2 3 ) = 1, a clear violation of 
Definition |2J 

We end this section remarking that Dcfinition|21is com- 
pletely defined only for finite chains. For infinite chains 
(N — > oo) one would have to calculate all G(n, 11,12) 
(and G(n,ii,i2, ■ ■ ■ ,i n -i)) to completely characterize a 
genuine 77-partite entangled state. Finally, the previ- 
ous definition does not imply that all genuine MES 
must have Tr {pf) < 1/2, Tr^-^J < 1/2, . . ., 

^ (Pjj+ii,— J+in-i) — ^/^ - ^ ^ S tml S 0n ^ y a sumc i en t 

condition for a state to be a genuine MES. 



B. Infinite Chains 

Currently there is an increasing interest on the rela- 
tion between entanglement and Quantum Phase Tran- 
sitions occurring in infinite spin chains |l6l ITtI I30L l3ll 
I35L l3rl Wq). For spin chains presenting a second order 
quantum phase transition (QPT) the correlation length 
goes to infinity at the critical point, thus suggesting in- 
teresting entanglement properties for the ground state 
of such models. Particularly interesting is the ID Ising 
model 33] , which is translationally invariant and presents 
a ferromagnetic-paramagnetic QPT. As we have seen in 
Sec. II, the generalized global entanglement is easily eval- 
uated for a system with translational symmetry. In this 
perspective, for the ID Ising model ground state, here 
we compute G(l), which is shown to behave similarly to 
the von Neumann entropy calculated in Ref. |30l |. and 
G(2, ii) for some values of i\. 

The ID Ising model with a transverse magnetic field is 
given by the Hamiltonian 



N N 

H = \J2^^ + i+z2' 



(37) 



This model has a symmetry under a global rotation of 
180° over the z axis (a x — > — a x ) which demands that 
(a x ) = 0. However as we decrease the magnetic field, in- 
creasing A, this symmetry is spontaneously broken (in the 



thermodynamic limit) and we can have a ferromagnetic 
phase with (a x ) 7^ 0. This phase transition occurs at the 
critical point A = A c = 1 where the gap vanishes and 
the correlation length goes to infinity. This transition is 
named quantum phase transition since it takes place at 
zero temperature and has many of the characteristics of 
a second order thermodynamic phase transition: phase 
transitions where the second derivative of the free energy 
diverges or is not continuous. It is worth noting that 
in the thermodynamic limit for A > 1 the ground state 
is two-fold degenerated. These two states have opposite 
magnetization. Here we will use the broken symmetric 
state for A > 1 and not a superposition of the two degen- 
erated states, which is also a ground state but unstable. 
For a more detailed discussion see Refs. [30ll3^. 

Now, let us explain how we can evaluate G(l) and 
G(2,ii) for the one dimensional Ising model. We need, 
then, the reduced density matrix of two spins, which is a 
4x4 matrix and can be written as 



Pij 



a,/3 



The coefficients are given by 



a/3 rp 



Pij 



(38) 



(39) 



and, as usual, Trjj is the partial trace over all degrees of 
freedom except the spins at sites i and j, erf is the Pauli 
matrix acting on the site i, a, (3 = 0, x, y, z where cr° is 
the identity matrix, and the coefficients p"? are real. 

Eq. (|39|l shows that all we need are the two-point spin 
correlation functions which, in principle, are at most 16. 
This number can be reduced using the symmetries of the 
Hamiltonian (|37|l . The translational symmetry implies 
that depends only on the distance \i — j\ = n between 
the spins so that we have pf? = p" 13 and p" 13 = p^ a . All 
these symmetries imply that the only non-zero correla- 
tion functions are: p" a , P 0x = p x0 = P x , P° z = P z0 = P z , 
and p xz =p%*. 

First, let us show the diagonal correlation functions 
and the magnetizations, which were already calculated 
in Ref. |33|. For periodic boundary conditions and an 
infinite chain we have: 



i+ni 



ff(-l) 

9(0) 



ff(-2) 
<?(-!) 



g(-n) 

g(-n+l) 



5(ra-2) 5(n-3) 

5(1) 5(0) • 
5(2) 5(1) ■ 

g(n) 5(77- 1) • 



(40) 



5(-77 + 2) 

5(-n + 3) 
5(1) 



z\2 



5 (7l) 5 (-77) . 



(41) 



(42) 
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be written as 



and 



(1 



= 5(0). 



0, 

A- 2 ) 1 / 



A < 1 
A > 1 



with 



and 



g(ri) = l (n) + XI [n+ 1) 



l(n) = - 

IT 



dk 



cos (kn) 



l + A 2 + 2Acos(fc)' 



(43) 



(44) 



(45) 



(46) 



2\2 



G(i) = i - ( P x y - ( P z ) 



(47) 



We are now left with the evaluation of p^ z = pf x . This 
calculation was made in Ref. [34| where the authors ob- 
tained the off-diagonal, time and temperature dependent, 
spin correlation functions. In the paramagnetic phase 
(A < 1) the ground state has the same symmetries of 
the Hamiltonian which leads to p% z = 0. For the ferro- 
magnetic phase (A > 1) an explicit evaluation leaves us 
with an expression in terms of intricate complex integrals 
which are not straightforward to compute. For this rea- 
son we will use bounds for this off-diagonal correlation 
function. 

We can obtain an upper and lower bound for this cor- 
relation function by imposing the positivity of the eigen- 
values of the reduced density operator py . For the Ising 
model these bounds result to be very tight as we can see 
in Fig. [7| and depend on n. In Ref. hfi some of the re- 
sults here discussed were presented using zero as a lower 
bound. It is worth mentioning that since both G(l) and 
G(2, «i) are decreasing functions of the square of the cor- 
relation functions, a lower (upper) bound for the latter 
implies an upper (lower) bound for the former. 



0.3 



0.1 



Upper 


Bound 


Lower 


Bound 




As we have already shown G(l) is the mean linear en- 
tropy of one spin which, due to translational symmetry, 
is equal to the linear entropy of any spin of the chain. A 
similar related analysis was done by Osborne and Nielsen 
|30| for the von Neumann entropy instead of the linear 
entropy. As well as G(l), see Fig. |UJ the von Neumann 
entropy is maximal at the critical point 30]. At that 
time Osborne and Nielsen did not give much importance 
to this result since they suspected that the von Neumann 
entropy of one spin with the rest of the chain does not 
measure genuine MES. However, for a translational sym- 
metric state it is a reasonable good indication of genuine 
ME as we have shown in previous sections. (We have 
explicitly studied the linear entropy but the same results 
apply to the von Neumann entropy. We have adopted 
the former mainly due to its simplicity and relation to 
the Meyer and Wallach global entanglement 

Analyzing Eq. I|47|l we can understand why G(l) is 
maximal at the critical point (A = 1). As we explain 
in what follows, it is (cr x ) the main responsible for this 
behavior of G(l). For A < 1 we have \a x ) = 0. After 
the critical point, however, (a x ) ^ 0. Moreover, for A > 
1 Eq. (|44l) tells us that (a x ) is a monotonic increasing 
function of A and that (a x ) — » 1 as A — > oo. Therefore, 
since (cr z ) is negligible for large values of A and (a x ) m 1 
(See Fig. fSJ we must have G(l) approaching zero after 
the critical point. 




0.5 



1.5 



Figure 8: (Color online) Magnetizations p x — (a x ) (black/da- 
shed line) and p z — (a z ) (red/solid line) as a function of A. 

We now analyze G(2,ii). Using the Ising model sym- 
metries Eq. (I24|l reads, 



Figure 7: (Color online) Bounds for obtained by impos- 
ing the positivity of the eigenvalues of the reduced density 
operator pij. 

Since we have all the correlation functions at hand we 
proceed with the calculations of G(l) and G(2,i 1 ). Re- 
membering that for the Ising model p v = Eq. (|20|) can 



G(2,n) = l--[2(^) 2 



„xz\2 



(Pnf + (Pl V f + (Pn) 2 } 



(48) 



With Eq. (|48() we can evaluate G(2, n) for any value of n. 
In Fig. IHlwe have plotted G(l) and the bonds for G(2, 1). 
We can see that both G(l) and G(2, 1) are maximum at 
the critical point A = 1. Notice that the bounds are 
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very tight and can barely be distinguished just in a small 
region for A > 1. Furthermore, G(2, 1) is always smaller 
than G(l), contrary to what was obtained using zero as 
a lower bound As well as in the case of G(l) we 

can see that the reason for G(2, 1) being maximal at the 
critical point is due to the behavior of (a x ) since it is the 
only function in Eq. I|48|) that does not change smoothly 
as we cross the critical point (see Fig. ^|for the other 
correlation functions). 
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Figure 9: (Color online) G(l) (red/dashed line) and the 
bounds for G(2,l) (black/solid lines). Note that they are 
maximum at the critical point. 
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Figure 10: (Color online) Two point correlation func- 
tions: pf x (red/solid), — p\ v (black/long-dashed), and p\ z 
(blue /short-dashed) . 

We have also plotted G(2,n) for n = 1, 7, and 15 
(Fig. Ill|) . We can observe that all of them are maxi- 
mum at the critical point and increase as a function of 
n (In Fig. ^] we have plotted only the upper bounds 
since the lower bounds produce very similar curves). We 
also note that G(2, 7) is very near G(2, 15) showing that 
G(2, n) rapidly saturates to a fixed value. At the critical 
point we have limn^oo G(2, n) = 0.675. This behavior for 
G(2,n) points in the direction of the existence of multi- 
partite entanglement at the critical point since any two 
spins are entangled with the rest of the chain and this 
entanglement increases with the distance between them. 
It is also interesting to confront this result with the fact 
that two spins that are separated by two or more sites 



are not entangled since their concurrences are zero |3l| . 
The behavior of the concurrence (G(n)) can also be un- 
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Figure 11: (Color online) G(2,n) for n = 1,7, and 15. From 
bottom to top n = 1,7, and 15 

derstood if we note that it can be expressed in terms 
of the one and two point correlation functions. While 
for the non-symmetric (ferromagnetic) state the analyti- 
cal expression for the concurrence is cumbersome for the 
symmetric one it is very simple. Fortunately, for the Ising 
model it was show that the concurrence does not change 
upon symmetry break |20l l2l| and it turns out to be 

c(n) = \(-i- p yy+ P T+Pn)- (49) 

From this expression we can see that the concurrence 
(Fig. Il2"|l does not depend on either the off-diagonal cor- 
relation function p^ z or on the one point correlation func- 
tions (magnetizations) . This is an interesting feature and 
helps us to understand why the concurrence is not max- 
imum at the critical point. 
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Figure 12: Concurrence for nearest neighbors. 



IV. CONCLUSION 

A TV-partite quantum system may be entangled in 
many distinct ways. To characterize and to define a 
good measure of entanglement for those systems is a hard 
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problem. The only simple alternative, valid whenever 
the joint iV-system state is pure, is to split the system 
into two partitions and compute the entanglement in that 
way. This bipartition could be constructed in many dif- 
ferent forms and thus give distinct amount of entangle- 
ment. One possible approach is to divide the system into 
two blocks of L and N — L subsystems and to compute 
the block entanglement [161 Il7| between the two blocks. 
However one could think of a situation where all of the 
subsystems in the block L are entangled with each other, 
as well as the subsystems of block N — L, but without any 
entanglement between the two blocks. For this situation 
the block entanglement would quantify a zero amount 
of entanglement, which is clearly not true. A valid bi- 
partition approach, which would be able to quantify the 
entanglement in such a situation, is to compute the en- 
tanglement for all kinds of bipartition and then to aver- 
age these to give the total amount of entanglement in the 
system. 

In this article we have formalized an operational multi- 
partite entanglement measure, the generalized global en- 
tanglement (Eq), firstly introduced in Ref. [Io|. For 

(n) 

n = 1 , E G recovers the Meyer and Wallach global en- 
tanglement measure 01 • However for n > 1 Eq^ to- 
gether with the auxiliary function G(n, i\, Z2, . . . , i n -i) 
quantify entanglement in the many distinct forms it is 
distributed in a multipartite system. We have shown that 
for some multipartite systems the original global entan- 
glement is not able to properly classify and identify mul- 
tipartite entanglement in a unequivocally way, whereas 
higher classes (n > 1) of E^ are. A genuine fc-partite en- 
tangled state is the one that cannot be written as a prod- 
uct \<p)i (g> j) of state vectors for any / < k, meaning 
that there is no other reduced pure state out of the joint 
fc-systems state. To completely quantify and classify the 
multipartite entanglement in this kind of state one would 

(n) 

have to compute all the E G classes up to n = k—1. How- 
ever we have observed that lower classes of E G , such as 
Eq^ and Eq^ , are sufficient to detect multipartite entan- 
glement. The computation of higher orders of Eq and 
of the auxiliary functions G(n,i\,i2, ■ ■ ■ ,in—i) is neces- 
sarily required only to distinguish and classify the ways 
the system is entangled. Although the calculation of all 
those higher orders may be operationally laborious it is 
straightforward to perform for finite N systems. Thus we 
have demonstrated for a variety of genuine multipartite 
entangled qubit states 0>HII that E G ^ and G(2,i\) are 



able to properly identify and distinguish them whereas 
Eq^ fails to do so. Inspired by the common characteris- 
tic presented by all G(2, i\) for those paradigmatic states 
we then discussed an operational definition of a genuine 
multipartite entangled state 0> E3 • 

Finite multipartite systems are interesting for funda- 
mental discussions on the definition of multipartite en- 
tanglement. Infinite systems on the other hand are inter- 
esting since multipartite entanglement may be relevant to 
improve our knowledge of quantum phase transition pro- 
cesses occurring in the thermodynamical limit. We have 
demonstrated that for the ID Ising model in a transverse 

(2) 

magnetic field both E G and G(2, i\) are maximal at the 
quantum critical point, suggesting thus a favorable pic- 
ture for the occurrence of a genuine multipartite entan- 
gled state. Moreover, the behavior of G(2, i\) and thus 

(2) 

E G can be easily understood as contributions of the one 
and two-point correlation functions giving us a physical 
picture for the behavior of the multipartite entanglement 
during the phase transition process. 

In conclusion the generalized global entanglement we 
presented has the following important features: (1) It is 
operationally easy to be computed, avoiding any mini- 
mization process over a set of quantum states; (2) It has 

(n) 

a clear physical meaning, being for each class E G the 
averaged n-partition purity; (3) It is able to order dis- 
tinct kinds of multipartite entangled states whereas other 
common measures fail to do so; (4) It is able to detect 
second order quantum phase transitions, being maximal 
at the critical point. (5) Finally, for two-level systems it 
is given in terms of correlation functions, and thus easily 
computed for a variety of available models. We hope that 
this measure may contribute for both the understanding 
of entanglement in multipartite systems and for the un- 
derstanding of the relevance of entanglement in quantum 
phase transitions. 
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